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Abstract The effects of dynamic localization in a solid-state system - a quantum dot - 
are considered. The theory of weak dynamic localization is developed for non- 
interacting electrons in a closed quantum dot under arbitrary time-dependent 
perturbation and its equivalence to the theory of weak Anderson localization is 
demonstrated. The dephasing due to inelastic electron scattering is shown to 
destroy the dynamic localization in a closed quantum dot leading to the classi- 
cal energy absorption at times much greater than the inelastic scattering time. 
Finally a realistic case of a dot weakly connected to leads is studied and it is 
shown that the dynamic localization may lead to a drastic change of the shape of 
the Coulomb blockade peak in the dc conductance vs the gate voltage. 
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1. Introduction 

The process of energy absorption by a quantum system with a time-dependent 
Hamiltonian underlies a large part of modern physics, both fundamental and 
applied. A generic Hamiltonian can be written in the form: 



H{t)=H Q + V q b{t), 



(1.1) 



where we explicitly separated the time-independent part Hq and the external 
perturbation V with the time dependence specified by a given function 4>{t). 

In the textbook example of the classical Drude absorption, qb{t) = E(t) is 
the time-dependent electric field which is often considered to be a harmonic 



2 



function of time E(t) = E w cos(ujt). Then the energy absorption rate Wq = 
H is given by the classical Joule heat formula: 

Wo = VE%<t{u), (1.2) 

where ct(uj) is the frequency-dependent conductivity and V is the volume of 
the system. In this classical picture the energy £ of a closed electron system 
grows linearly with time. Assuming the model of non-interacting electrons 
with the density of states v related by vV = 1/5 with the mean separation 5 
between discrete one-particle levels we find: 

/ds 
e[/00 - M £ )\ y«y + const, (1.3) 

where /(e) is the (non-equilibrium) electron energy distribution function, M £ ) = 
9{—e) is the zero-temperature Fermi-step, and T e g is the effective temperature 
of electron system. One concludes from Eq.(1.3) that the classical absorption 
picture corresponds to diffusion in the energy space: the width of the electron 
energy distribution T e g increases with time according to the diffusion law: 

T? s = D t, (1.4) 

where Do = Wq8 is the energy diffusion coefficient. 

Apparently, this corresponds to the Markovian random walk over the spec- 
trum of non-perturbed system caused by absorption and emission of energy 
quanta Huj, each absorption and emission events being independent of other 
ones. Note that in this random walk picture the only limitation on the valid- 
ity of Eqs.(1.2,1.4) is the finite width of the quasi-continuous energy band that 
leads to a saturation in the effective temperature or absorbed energy. The trivial 
example of such a saturation are Rabi oscillations in a two-level system where 
the averaged energy does not change with time. 

In the past two decades attention of the scientific community was drawn to 
a different and much less trivial example of saturation in the time-dependent 
energy of closed driven systems called dynamic localization (DL). In this case 
the spectrum of Ho is essentially unlimited in the energy space, yet, after a 
certain time the absorption rate W(t) = D(t)/5 or the energy diffusion coef- 
ficient D(t) in Eq.(1.4) vanishes. The DL in the energy space was observed 
in numerical simulations on the kicked quantum rotor (KQR) - particle on a 
circle with Hq = —d 2 /d9 2 and <j)(t) being aperiodic sequence of (5-pulses [1], 
as well as in an actual experimental realization of the KQR - trapped ultracold 
atoms in the field of a modulated laser standing wave [2]. The mapping of 
the KQR to the quasi -random Id Anderson model has been done in Ref. [3], 
a similar analogy was exploited in Ref. [4] to demonstrate the DL in a meso- 
scopic disordered ring threaded by a magnetic flux growing linearly in time. In 
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Ref. [5] an analogy between the KQR and band random matrices was pointed 
out, the latter have been reduced to a Id nonlinear a model [6]. In Ref. [7] 
the direct correspondence between the KQR and a Id nonlinear a model was 
demonstrated. 

The physical origin of dynamic localization in the energy space of driven 
systems is essentially the same as for its real space counterpart - the Anderson 
localization - in stationary systems. This is the quantum nature of absorption 
and emission which does not reduce merely to a quantized step hw of the ran- 
dom walk. According to basic principles of quantum mechanics one should 
consider paths in the energy space between an initial and a final point each 
of them consisting of many such steps. The transition probability is given by 
a square modulus of the sum of corresponding amplitudes where all interfer- 
ence terms have to be taken into account. It is these interference terms that 
makes the picture of independent absorption and emission events incomplete 
and eventually may lead to dynamic localization. In particular the weak dy- 
namic localization (WDL) exhibits itself as a time-dependent negative correc- 
tion 5D(t) to the energy diffusion coefficient which magnitude is controlled by 
a small parameter S/T <C 1, where V = (V 2 }/5 is the typical radiation width 
of energy levels. 

The present study consists of three parts. In the first part we develop an 
analytical theory of weak dynamic localization in a closed system of non- 
interacting electrons driven by a time-dependent perturbation with an arbi- 
trary function 4>(t). To accomplish this goal we use a variant of the diffuson- 
cooperon diagrammatic technique [8, 9] in the time domain that has been de- 
veloped on the basis of the Keldysh formalism for non-equilibrium systems. 
We focus on the case - most relevant for the quantum dot application - where 
both the unperturbed system Hq and the perturbation operator V are described 
by the Gaussian ensemble of random matrices or by an equivalent random field 
Hamiltonian considered in the zero-mode approximation. Such description has 
proved [see [8] and references therein] to be valid for low-energy domain of 
diffusive quantum dots and is believed to be also accurate for ballistic quan- 
tum dots with irregular boundary. We establish an intimate relation between 
WDL and dephasing by time-dependent perturbation using the notion of no- 
dephasing points first introduced in Ref. [10]. We show that the existence and 
character of weak dynamic localization depends crucially on the spectral and 
symmetry properties of 0(i). Furthermore we show that by merely changing 
(j)(t) one can obtain the dependence of SD(t) that reproduces the dependence 
of WL correction to electric conductivity 8a {tip) on the dephasing time in 
quasi-one dimensional wires, two- and three-dimensional systems. 

In the second part we give a brief description of the effect of electron- 
electron interaction on dynamic localization in closed systems assuming that 
in the absence of interaction the strong dynamic localization occurs at t > i* 
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with the effective temperature at saturation T* = \fD^tl. We show that in 
contrast to the Anderson localization where the hopping conductance is finite 
but exponentially small at low temperatures, the strong dynamic localization is 
destroyed by electron-electron interaction even in the case where the character- 
istic time of inelastic electron-electron scattering t ee = t ee (T*) is much larger 
than the localization time t*. What is left of DL in closed systems with electron 
interaction is a suppression of the energy absorption rate W(t)/Wo ~ t^/t cc 
at intermediate times, t* < t < t ee . For times longer than t ee ln(t ee /t*) the 
classical absorption rate Wq is again recovered. 

Finally, in the third part we address a realistic case of a quantum dot in a 
steady-state regime under time-dependent perturbation where both electron- 
electron interaction and electron escape into leads are taken into account. We 
show that a signature of dynamic localization can be observed in an almost 
closed quantum dot with the escape rate 7 esc <C S. This is a plateau at the tail 
of the Coulomb blockade peak in the dc conductance vs. the gate voltage. We 
obtain an analytic expression for the shape of the Coulomb blockade peak in 
driven quantum dots and identify a region of parameters where the plateau can 
occur. The main results of the first two parts are published in Ref.[ll, 12] and 
a detailed study of the Coulomb blockade regime in a driven quantum dot is 
presented in [22]. 

2. Weak dynamic localization 

The goal of this section is to show how the cooperon-diffuson diagrammatic 
technique [14] which is the main theoretical tool for describing weak Anderson 
localization (WAL) and mesoscopic phenomena, can be extended as to include 
non-equilibrium processes considered in the time-domain. This consideration 
is based on the Keldysh technique [13] and is described in detail Ref.[8, 9]. 
Here we present only the main practical hints which allow to derive the WDL 
correction along the same lines as those used to derive WAL corrections to 
conductivity [14]. 

We start by the standard [13] expression for the energy distribution function 

/(£, t) = \ ~ \ J h (t + I * - I) dr, (1.5) 

in terms of the Keldysh Green's function Q K (t, if): 

g K (t,t') = -2irivh(t,t'). (1.6) 

Then the time-dependent absorption rate W(t) = D(t)/5 which is related to 
the time-dependent energy diffusion coefficient D(t), is given by: 
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Next we note that the Keldysh Green's function for the free electron gas cou- 
pled with an external time-dependent field by H e f(t) = V <j)(t) can be ex- 
pressed in terms of retarded Q R and advanced Q A Green's functions as follows: 

G K (t, = j dhdt', g R (t, h) g A (t' 1 ,t') ho(h - ti) [Kef (ti) - 7U(h)] 
+ J dt 1 h (t-t 1 )g R (t 1 ,t')- J dt 1 g A (t,t 1 )h (t 1 -i?),(i.«) 

where ho(t) is the Fourier-transform of the equilibrium distribution function 
tanh(e/2T) = 1 — 2fo(e) which is supposed to hold in the absence of per- 
turbation 7i e {. Note also that the functions g R ( A \t,t') take account of the 
external time-dependent field in all orders in Tt e f and therefore depend on two 
time arguments and not only on their difference. 

The perturbation operator in Eq.(1.8) is not necessarily random and without 
loss of generality it can be assumed to be traceless TrH e { = 0. For instance 
one can consider the perturbation by a time-dependent space-homogeneous 
electric field E(i) = — d t A = E o 0(i), where Ti c { = — evA(i) with zero 
average of the electron velocity v over the Fermi surface. 

2.1 Diffuson-Cooperon diagrammatic technique in the 
time domain 

Now we assume that the Hamiltonian Hq of the unperturbed system is that 
of free electrons in a Gaussian random impurity field U(r) which results in 
the diffusive electron dynamics with a small elastic mean free path £ <C L 
compared to the system size L. Averaging over disorder realizations U(r) 
can be done using impurity diagrammatic technique [15]. For relatively weak 
external fields such that 

eE £ <C huJ (1.9) 

the essentially nonlinear effect of the field is on the two-particle correlation 
functions -diffusons and Cooperons- while for the averaged single-particle 
Green's function an expansion in H e { up to the second order is sufficient: 



{g R{A) (t, 0) = s(t - 1 1 ) [g r ^ + g r ^ w rf (t) g r ^ 

+ G R ^ H ei (t) G R ^ Kef (t) G R ^ 



(1.10) 



In Eq.(l.lO) it is assumed that the elastic scattering time tq = £/vf is the 
smallest relevant time scale, and 
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where £ p = ^ - e F . 

Using Eq.(l.lO) one can decouple the time- and momentum- integrations. 
The latter for a closed loop of G R ^ functions can be done as a simple pole 
integral v J d£ G R ^ (£)...G R ( A ) (£), as in the stationary case. In particular, 
one can repeat the standard derivation of the diffusons and Cooperons as the 
corresponding ladder series of the impurity diagrammatic technique. 

The diffuson 

(g R {v,v';t + ,t' + )g A (r',v-t'_,t-)) (1.12) 
= 2ttu 6{t - t') e irKd{t ' T) D T (t, t'; r, r ') e -^K d (t',r) ) 

where t± = t±r/2, t'± =t'± r'/2, K d(c) (t, r) = A(t + r/2) T A(t - r/2), 
obeys the equation [8]: 



d o 



xD T (t, t'; r, r) = 5(t - t')S(r - r') 



(1.13) 



where D = vpijd is the diffusion coefficient in the d-dimensional space. 
The corresponding equation for the Cooperon 

(g R (r,r';t + ,t' + )g A (r,r';t-,t'_)) (1.14) 
= 2nu 6{t - t') e irKc( '' r) C t (r, r'; r, Y >y^^') 



reads [8]: 



d o 
2— - DV 2 r - irE 
or 



(*+§)-*(«-§) 



xC t {T,T' ;r,r') = 2S(t - t')S{t - r') (1.15) 

Equations (1.13,1.15) should be supplemented with the Neumann boundary 
conditions: 

dD T (t,t';r,r') dC t (r,r';r,r') 

a = U ' 5 = U - (1-16) 

an an 

In contrast to the initial transverse gauge, in the longitudinal gauge we switched 
to in Eqs.(1.12-1.15) the boundary conditions do not depend on the time- 
dependent perturbation. 

Without time-dependent perturbation Eqs. (1.13,1. 15) are just simple diffu- 
sion equations which have a complete set of stationary solutions <^(r) such 
that DV^Lp^r) = E^ip^r) including the zero mode i^o = const. One can 
use this set as a convenient basis for solving Eqs.(1.13,1.15) in the presence of 
time-dependent perturbation [8]. In the limit 



eE L, hio < E Th , 



(1.17) 
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where Eth = hE\ ~ hD/L 2 is the Thouless energy, the zero mode makes 
the main contribution, other modes should be treated as perturbations. Then 
the second-order perturbation theory yields for the zero-mode diffuson and 
Cooperon [8, 9, 16]: 



D T (t, = 0(t - t>) exp j-r J* [<f> (t" + I) - (V - I)] 2 dt'j 

(1 

C t (r, r') = 0(t - r') exp |-| jT (t + \) - (f - \)] ' di/} , 



(1.18) 



where 



= [erp^Ep] 2 ^ (e£ L) 2 



(1.19) 
(1.20) 



and r 0M = fdrfo r<^(r). 

The validity condition for the zero-mode approximation can be obtained 
using Eqs.(1.9,l. 17,1.20): 



£ 
L 



< 1. 



(1.21) 



It is remarkable that the same expressions for the diffuson and the Cooperon 
can be obtained [8, 9] if one starts from the Hamiltonian Eq.(l.l) with Ho and 
V taken from the Gaussian Orthogonal ensemble of real symmetric N x N 
random matrices with the probability distributions: 



Vh °c exp 



4N5 2 



Vy oc exp 



irtr V 2 ' 
4T5 



(1.22) 



This demonstrates that for small conductors (quantum dots) and relatively 
weak external fields where Eq.(1.21) satisfies, the non-random time-dependent 
perturbation Vo4>(t) = erEo 4>(t) in a random system is equivalent to a random- 
matrix perturbation V with (V 2 ) ~ V 2 '/ 1 g, where g = E^k/d ^ 1- This 
statement is believed to hold not only for diffusive quantum dots but also for 
generic ballistic quantum dots. 

However, perhaps the simplest way to generate a correct diffuson-cooperon 
diagrammatic technique in the zero-dimensional limit Eq.(1.21) is to consider 
in Eq.(l.l) 

H = ^- + U(r) + V(r)<f>(t) (1.23) 
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where U(r) and V(r) are independent Gaussian random fields with zero mean 
values and the following correlation functions: 

(U(r) U(r')) = — L_ 5(r - r'), (V(r)V(r')) = —S(r - r'); Tr « 1. 
ZTTUTq iris 

(1.24) 

We will use this representation below because it corresponds to a minimal 
deformation of the impure-metal Hamiltonian that is a starting point of a con- 
ventional impurity diagrammatic technique and because it helps to avoid subtle 
effects of the boundary conditions Eq.(1.16) which are important for a deter- 
ministic global perturbation. 

2.2 One loop DL correction 

As a first step of averaging over disorder in Eq.(1.8) we single out the ladder 
series that represent diffusons and Cooperons. As usual [14] the diagrams are 
classified with respect to the number of diffuson or Cooperon loops. The zero- 
and one-loop diagrams are shown in Fig. la and Fig.lc, respectively. The next 
step is to properly average the remaining products of retarded and advanced 
single-particle Green's functions (denoted by R and A in Fig.l). So we obtain 
the Hikami boxes denoted by the shadowed polygons in Fig. lb and Fig. Id. In 
order to compute the Hikami boxes, one substitutes Eq.(l.lO) for the averaged 
retarded or advanced Green's functions in Fig.la,c retaining only terms second 
order in the perturbation operator H e { = V(r) 4>(t). Then one performs aver- 
aging over the random field V(r). It is shown in Fig.2 as the dashed-dotted line 
connecting two vertices. In addition to that rules of the impurity diagrammatic 
technique require additional averaging over the impurity field U(r) which is 
shown by the dotted line in Fig.2. The further calculation of the Hikami boxes 
is standard and reduces to pole integrals of G R ( A \£) functions Eq.(l.ll). 

One can easily see that the diagrams Fig.2a and Fig.2b have opposite signs. 
As the result the Hikami box Fig.2a,b in the zero-loop diagram Fig. lb is pro- 
portional to [4>(t'+T/2)-<p(t'-T/2)] 2 ,so that the function h(t+T /2,t-T /2) 
in Eq.(1.6) obeys in the zero-loop approximation the following equation: 

Wt+I,t-I)=0. (1.25) 

If one is interested in the energy resolution Ae » u, where a; is a typical 
period of a time-dependent field, it suffices to consider lot <C 1. Then from 
Eq.(1.25,1.5) we obtain diffusion in the energy space: 

- r (3, 4>(f)?^} /(£,() = 0. (1.26) 



"5 W'-5. 
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D T (t,t') 



b) 



t'+T/2 tH~ 




\ 1+ T/2 



C) 



D T (t,t') 



-o 



h(x)^) 



C t ,^ /4 (^/2,-^/2) 



d) 



Figure 1.1. Diagrams for the Keldysh Green's function Q K (t+r/2, t — r/2): a), the zero-loop 
diagram describing the classical absorption; the diffuson is given by a ladder series, the retarded 
or advanced Green's functions are depicted by solid lines b). the same diagram with Hikami 
boxes shown by shadowed loop and a diffuson denoted by a wavy line; c). one loop WDL 
correction with a diffuson and a Cooperon given by ladder series; d). one loop WDL correction 
with the Hikami box denoted by a shadowed rhomb; the diffuson and the loop Cooperon are 
denoted by the wavy lines 



If in addition we are interested in the time resolution At » uj 1 , the time- 
independent diffusion coefficient can be defined: 



d = T{d t <p{t)) 2 -rv 2 , 



(1.27) 
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Figure 1.2. Hikami boxes: solid lines denote retarded or advanced Green's functions, 
Eq.(l.ll); the dashed-dotted line and the dotted line represent the (V V) and the (U U) cor- 
relator Eq.(1.24), respectively. The black circle denotes the V<f>(i) vertex in Eq.(1.8); the open 
circle corresponds to the V<f>(t) vertex in Eq.(l.lO). 

where (dt4>) 2 averaged over time interval much larger than the typical period 
of oscillation is supposed to be well defined and independent of time. 

One can see from Eq.(1.26) that the parameter T defined in Eqs.(l. 20,1.22,1. 24) 
has a physical meaning of an inverse time of making a step ~ u in a random 
walk over the energy space. 

The correction 5D(t) to the energy diffusion coefficient is described by the 
diagrams Fig.lc,d where the Hikami box is given by diagrams Fig.2c-g. One 
can easily see that in the zero-mode approximation (zero external momentum 
in any vertex of the rhomb) three diagrams Fig.2c-e cancel each other. The 
same is valid for the diagrams Fig.2c'-e'. The remaining diagrams Fig.2f,g 
have opposite signs, so that the entire Hikami box in Fig. Id is proportional to 
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(<f)(t' + r/2) - <p{t' - r/2)) {<p{t' - £) - 4>(t' - £ - r)). Then Eq.(1.7) yields 
[11]: 

= — / d t 4>(t) d t <p(t - o c t _^ -0 de, (1.28) 

^ JO 

where C t _£/ 2 (£> ~C) is tne zero-mode Cooperon given by Eq.(1.19) and t is 
time passed since the onset of the time-dependent perturbation cf>(t). Note that 
the correction SD(t) is a purely quantum mesoscopic effect that vanishes if 
the mean level spacing 5 between energy levels of the non-perturbed system is 
equal to zero. 

Equation (1.28) that describes the weak dynamic localization is the main 
result of this section. Below we establish some of its implications. 

2.3 No-dephasing points and dynamic localization 

Let us consider the simplest example of a harmonic perturbation <p{t) = 
cos ut. Then for u£, S> 1 one can write: 

C t - € /2(£, -0 = exp j-2£r sin 2 Lt - ^fj \ . (1.29) 

A remarkable phenomenon is that the negative exponential that describes de- 
phasing by a harmonic perturbation vanishes on a certain set of no-dephasing 
points: 

u}£ k = 2oot - 2irk, k = 0,±1,±2, ... (1.30) 

At large T£ only a close vicinity of these points contributes to the integral 
Eq.(1.28). Expanding sin 2 around no-dephasing points and performing the 
Gaussian integration we arrive at: 



where 




(1.31) 



(1.32) 



This result is valid at t <C i* and £ ~ t S> 1/T which is only compatible if 

T > 5. (1.33) 

Remarkably, the WDL correction Eq.(1.31) increases with time. Equation 
(1.31) sets a time scale t* where the negative WDL correction is of the or- 
der of the classical energy diffusion coefficient and one may expect the strong 
dynamic localization. 
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Note that it is only because of the no-dephasing points which are zeros of 
the dephasing function 



(1.34) 



that the correction 5D(t)/Do grows with time and may become of order one 
despite a small parameter 5/F. For instance, the white-noise perturbation 4>(t) 
such that the time average 4>(t)(f)(t') oc S(t — t'), results in a constant dephasing 
function T c = T^. Then the WDL correction Eq.(1.28) is finite SD(t) ~ 
Dq 5/T at t —> oo and small compared to Dq. 

In order to understand the physical meaning of no-dephasing condition we 
consider two electron trajectories in real space with loops traversed in opposite 
directions (see Fig.3), the traversing time being t. Interference between such 
trajectories is known to be a cause of the weak Anderson localization. Let us 
assume that a time-dependent vector-potential A(t) is present. Then the phase 
difference between trajectories 1 and 2 in Fig.3 becomes random: 

5<p = f A(t') MO - v 2 (t')\ dt'. (1.35) 
Jo 

Using the condition vi(t') = —V2(t — if) we obtain 

ft/2 

5ip= / [A(t/2 + t') + A(t/2-t')]v 1 (t' + t/2)dt'. (1.36) 

J-t/2 

If the typical period of oscillations in A(t) is large compared to the velocity 
correlation time, one obtains after averaging over directions of velocity v\ : 

rt/2 

((^) 2 )oc / [A(t/2 + t') + A(t/2 -t')] 2 dt'. (1.37) 

J-t/2 

For a harmonic A(t) = smut it is easily seen that the phase difference is zero 
for trajectories of any shape provided that the traversing time t is synchronized 
with the period of the field: 

2irk 

t = , k = 1,2, .. (1.38) 

that is, the dephasing is absent for trajectories which traversing time is equal 
to a multiple of the period of perturbation. 

2.4 No-dephasing points and generalized time-reversal 
symmetry 

Now let us consider a generic periodic perturbation: 

oo 

4>(t) = ^2 A n cos(mujt + (p m ). (1-39) 

m=l 
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Figure 1.3. Time-reversal conjugate trajectories. 
First of all we require the diffusion coefficient Eq.(1.27) to be finite: 

p oo 

D = - J2 m2A2 m<™- (i" 40 ) 
m=l 

This is only so when decreases faster than 1/m 3 . In particular, it is not 
the case when <fi(t) is the periodic (5-function as in the KQR model [1, 3]. In 
similar cases the sum Eq.(l-40) should be cut off at m ~ E^h/uj, so that 
the energy diffusion coefficient Dq grows with increasing Exh or the size of 
matrix N ~ E^hjb in the equivalent random matrix model Eq.(1.22). 

The cooperon dephasing function Eq.(1.34) that corresponds to Eq.(1.39) 
takes the form: 

oo 

r c (t) = r A l sm 2 (W + <Pm)- (1.41) 

m=l 

The no-dephasing points tk determined by F c (tk) = exist only provided 
that phases of all harmonics are sinchronized: ip m = imp (see Fig.4a). This 
condition implies that the perturbation (f>(t) can be made a symmetric function 
by a shift of the origin <j>(t) = 4>(t — <p/uj): 



(1.42) 




Figure 1.4. The Cooperon dephasing function T c (t) for (a) periodic <f>(t) obeying Eq. (1.42): 
a regular array of zeros; (b) generic periodic <j>(t): a gap; (c) quasi-periodic <j>{t) with two 
incommensurate frequencies: a pseudo-gap. 



Apparently, this is a generalization of the time-reversal condition for the case 
of time-dependent perturbation and for quantities that survive averaging over 
time intervals much larger than the period of oscillations in 4>(t). It has been 
first established in Ref.[10]. 

If this condition is violated (see Fig.4b) and minr c (t) = 7 > 0, the 
Cooperon Cf(£, — £) acquires an exponentially decreasing factor e™ 27 ^ that 
limits the growth of WDL correction for tj S> 1. For a generic case 7 ~ T 
we have SD(t)/Do ~ {S/T) *C 1 and the one-loop WDL correction given by 
Eq.(1.28) can be neglected. 

The analogous situation arises in the weak Anderson localization when the 
time-reversal symmetry is violated by a constant magnetic field or magnetic 
impurities. In this case the one-loop WAL correction also gradually vanishes 
leading to the anomalous magneto-resistance. The difference is that in the case 
of WDL the Hamiltonian remains real, represented e.g. by the orthogonal 
ensemble of random matrices. Yet the time-reversal symmetry can be broken 
just by time-dependence of perturbation. From this perspective the harmonic 
perturbation is an important exception rather than a paradigm of a periodic 
perturbation. 
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The analogy with the WAL suggests that the main WDL correction in the 
case where the condition Eq.(1.42) is violated is given by the two-loop dia- 
grams containing only diffusons. Indeed, let us consider the diffuson dephas- 
ing function: 



r d (r)=r[0 (f + I)-^'-!)"' 2 



2 

For cf)(t) given by Eq.(1.39) it reduces to 



(1-43) 



oo 

r d (t) = 2r^^sin 2 (m^). (1.44) 

m=l 

As was expected, Td(t) is independent of phases (p rn and has a set of no- 
dephasing points = 2nk/oj for any periodic perturbation Eq.(1.39). As 
the result, the two-loop WDL correction is a growing function of time [17]: 

^M = _*A, ( , 45) 

D Q 24 U ' 

where 

U = -^p- (1-46) 

2.5 Incommensurability of frequencies, avoided 

no-dephasing points and weak Anderson localization 
in higher dimensions. 

Now let us consider the perturbation <f>(t) as a sum of several incommensu- 
rate harmonics: 

d 

<p(t) = ^2 A m cos(uj m t + <p m ). (1.47) 

m=l 

In this case the Cooperon dephasing function T c (t) is given by: 

d 

r c (t) = r A l sm 2 (Lo m t + <p m ), (1.48) 

m=l 

where all uj, m do not have a common multiplier. In this case one has a distribu- 
tion of avoided no-dephasing points as in Fig.4c. 

In order to compute the one-loop WDL correction Eq.(1.28) we expand each 

of the exponentials [sf. Eq.(1.29)] exp j-2£rM^ sin 2 (u m t - + ip m ^ } 
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in a Fourier series using: 

+00 

e z cos(2ui m t+2tp m ) _ J n (2) e in m (2uj m t+2tp m ) (149) 



n m =— 00 



where I nvn (z) is the Bessel function. 
Then Eq.(1.28) reduces to: 



m=— 00 tij=-oo m=l 

x JJ / nm (z m )e _2m e inm(2Wm *- a;m$+2¥ ' m) , (1.50) 



where z m = £rA^ and the prime denotes the derivative d/dz m . 

In order to obtain the behavior of 5D(t) averaged over time intervals much 
larger than the typical period of oscillations one should take into account only 
n m obeying the constraint: 

d 

^2n m uj m = 0. (1.51) 

m=l 

It is exactly the condition of incommensurability of frequencies oj m that the 
constraint Eq.(1.51) can be satisfied only if all n m = 0. One can see that at 
such a condition 8D{t) does not depend on the phases ip m and is given by: 

■px rt d d 
8D(t) = -— ^ E u2 m A l ^h(z m ))' J] I (z m ) e~*™ (1.52) 

Further simplification is possible if we consider A m = 1 and Ft S> 1. In this 
limit (lnIo(z m )) = 1, Io(z m ) « e Zm / yj2irz m and one arrives at [11]: 

SD(t)_ 6 r dz 



D ttTJ^ (2irz) d / 2 ' 

Remarkably, the dependence on t in Eq.(1.53) coincides with the dependence 
of WAL corrections to conductivity of a d-dimensional conductor on the de- 
phasing time Ttp. In the case of two incommensurate harmonics the WDL 
correction grows with time logarithmically: 

For three and more incommensurate harmonics SD(t)/ 'Dq ~ S/T is small. It 
is yet to be studied whether or not a dynamic localization transition is possible 
for d > 3 at some r cr i t ~ 5 as it is claimed to be the case for the KQR [18]. 



Dephasing and dynamic localization in quantum dots. 
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3. Role of electron interaction on dynamic localization in 
closed quantum dots. 

So far we considered energy absorption in a closed system of non-interacting 
fermions under time-dependent perturbation. Although the many -body charac- 
ter of the system was important for making the diffusion in the energy space a 
mechanism of energy absorption due to the presence of the filled Fermi see, in- 
elastic processes of interaction were completely disregarded. We have demon- 
strated how the weak dynamic localization arises in such systems and how 
much does it resemble the Anderson localization (AL) corrections to conduc- 
tivity. This analogy makes us believe that WDL corrections for periodic per- 
turbations evolve with time to produce the strong DL in the same way as WAL 
corrections lead to a vanishing conductance of low-dimensional systems of 
sufficiently large size in the absence of dephasing mechanisms. So we assume 
that in a quantum dot without electron interaction at times t S> t* the energy 
diffusion coefficient becomes exponentially small. By this time the electron 
temperature rises from a low fridge temperature to the effective temperature 
T* - cjT/6. 

The question is what happens if electron interaction is present. More pre- 
cisely, what happens if the dephasing rate 7^ = 1 / t cc related with interaction 
is small but finite: 

< 7^*) « 1. (1.55) 

In the problem of strong AL the conductance a remains exponentially small 
provided that the phase-breaking length is much greater than the localization 
length 3> L\ oc and the temperature is much smaller than the mean level 
spacing in the localization volume 5\ oc = (i/Lf^ The first condition en- 
sures localization of eigenstates. The second one makes hopping over pre- 
formed localized states an exponentially rare event, so that a ~ e~~^ Sloc / T ^ a 
with some positive exponent a. 

The point of a crucial difference between the strong AL and the strong DL 
is that in the latter it is the inverse localization time i" 1 ~ 5 2 /V that plays a 
role of 5\ oc . In the case when S is the smallest energy scale we have 

1 

r* ~ -=- » -. (1.56) 

t* 

As a result, no exponentially small factor arises in the hopping over the Floquet 
states in contrast to the hopping over localized states in a real space. 

The qualitative picture that follows from the consideration [12] based on the 
Fermi Golden Rule is the following. A given electron absorbs energy from the 
external field during the time i* until the strong DL develops itself. At t > t* 
the dynamically localized state with the energy T c g ~ T* above the Fermi sea 
is formed which is characterized by the exponentially small absorption rate. 
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So it continues until after a time t cc a collision with another electron breaks 
the phase-tuning necessary for DL. The "efficiency" of each collision in the 
destruction of DL is close to 100% because the typical energy transfer is large 
T* 3> Then the absorption rate jumps to a classical value and it takes 

another t* to reach DL for the second time, but now with energy 2T*. So, for 
a given electron relatively short time intervals ~ of absorption are followed 
by larger periods ~ t cc of DL, the energy increasing by T* after each circle of 
absorption. The averaged absorption rate W- m is given by: 



W 



~\~ 



(1.57) 



Were electron-electron collision time t cc energy independent, the averaged ab- 
sorption rate would stay constant at all times t > t*. However, l/t cc typically 
increases with increasing the energy T c g, so does also the averaged absorption 
rate. Then we obtain from Eq(1.3): 



dt 



Dn 



U + t ee (T eS ) 



(1.58) 



The time t cc that corresponds to inelastic collisions of electrons in a quantum 
dot with the energy transfer ~ T c g- <C Exh has been calculated in Ref.[19]: 



(1.59) 



For t* <C t ee (T e s) Eqs.(1.58, 1.59) suggest that the effective electron temper- 
ature increases exponentially: 



T eff = T* exp 



2tee(-^*) 



(1.60) 



rather than remaining a constant ~ as it does for non-interacting electrons. 
The corresponding absorption rate Wi n (t) takes the form: 



Wi n (i) = W 



tee(T* 



exp 



t > U 



(1.61) 



The exponential growth of effective temperature continues until tee(T e g) de- 
creases down to t*. At larger times t > t^T*) ln(t ee (T*) /t*) the classical 
absorption rate is restored: W- m = Wo- The sketch of the time-dependence 
of absorption rate is shown in Fig.5. The main conclusion of the qualitative 
picture of Ref.[12] described above is that the dynamic localization is a tran- 
sient phenomenon that is destroyed at long enough times by an arbitrary small 
electron interaction. 
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Figure 1.5. Absorption rate vs time in a closed system with electron-electron interaction. 

Note, however, that the consideration of Ref.[12] is based on the Fermi 
Golden Rule. The applicability of this description is known [20] to be re- 
stricted in quantum dots by the region of relatively large electron energies 
T e fj 3> \JExh5- At smaller energies the localization in the Fock space should 
take place with the effect of vanishing l/t cc . Therefore it is likely that at 
T* ~ ujY/8 < \/Eth6 the dynamic localization survives electron-electron 
interaction. 



4. Dynamic localization in an open quantum dot and the 
shape of the Coulomb blockade peak. 

The transient character of DL and difficulty to realize a closed system in 
condensed matter physics because of the phonon heat transport leaves only 
one option for possible observation of DL: the steady state in an open system 
under periodic pumping. In an open system Eq.(1.57) should be generalized to 
include the dephasing due to electron escape into leads along with dephasing 
by inelastic electron collisions considered in the previous section: 

Wi n = Wo - , (1.62) 

where the energy-dependent dephasing rate is given by: 

7 V = 7csc + - L y (1-63) 

From Eq.(1.62) it follows that in order to observe any deviation from the clas- 
sical absorption rate Wo the escape rate 7 CSC must be much smaller than the 
inverse localization time 1/t* ~ S 2 /F. Then we obtain the necessary condi- 
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tion for DL: 




(1.64) 



The last inequality Eq.(1.33) is the basic approximation adopted in this study 



Eq.(1.64) shows that the DL effects considered could be observed only pro- 
vided that the quantum dot is almost closed, i.e. 2ip <^ i This is exactly 
the condition where the Coulomb blockade effects [21] are relevant. Under 
this condition the low-temperature linear dc conductance of the dot is strongly 
peaked near the value of the gate voltage V g = V c where the states with N 
and N + I electrons on the dot are nearly degenerate. At equilibrium condi- 
tions where the electron energy distribution inside the dot is Fermi-Dirac with 
the same temperature as in the leads, the width of the peak is proportional to 
the temperature and thus can be used as a thermometer. In the case where 
the dot is subject to time-dependent perturbation, the electron energy distri- 
bution inside the dot can be drastically different from that in the leads and 
a non-trivial question arises which of these two distributions affect the linear 
dc conductance. The answer for an open dot where the Coulomb blockade 
effects play no role is given in Refs.[8, 9]. It appears that in this case the lin- 
ear dc conductance is sensitive only to the electron energy distribution in the 
leads and thus cannot be used to probe energy dynamics in the dot. However, 
with the Coulomb interaction inside the dot taken into account and for the case 
T e s S> 5 the linear dc conductance becomes sensitive to the effective electron 
temperature T e g inside the dot. Below we consider only the extreme case of 
the Coulomb blockade. We also assume that the inelastic electron interaction 
is sufficiently strong to produce the Fermi-Dirac form of the electron energy 
distribution inside the dot, though with the temperature T e g much larger than 
that in the leads. We concentrate on the qualitative changes in the Coulomb 
blockade peak shape characteristic of the dynamic localization in the approxi- 
mation of sequential tunneling [21]. More detailed consideration including the 
role of inelastic co-tunneling and cooling by phonons is presented in Ref.[22]. 

The steady state of a dot under ac perturbation is determined by the global 
energy balance W- m = W out , where the pumping of energy by the ac per- 
turbation is governed by Eq.(1.62) and W out is the cooling rate. The most 
interesting case corresponds to a situation where the main cooling mechanism 
is due to electron escape into cold leads. Denoting by U = V g — V c the de- 
tuning of the gate voltage from the peak center and introducing the parameter 
x = U I (2T e fj) we obtain [22] for the case of leads maintained at a temperature 



of DL. 



T < T eff : 




(1.65) 



Dephasing and dynamic localization in quantum dots. 2 1 

where 7 = 2^f csc (U = 0). In the same approximation we obtain for the linear 




Figure 1.6. The function (G/Go) in Eq.(1.67). 

dc conductance of the dot: 

G(U) x tanhx 

Gq fn(2 coshx) 



(1.66) 



Remarkably, the r.h.s. of Eqs. (1.65, 1.66) depends only on the single variable 
x, so that one can establish a relationship similar to the Wiedemann-Franz law: 

W ont = (Jj T^W(G/G ). (1.67) 

Now, let us assume that 

7<A < 1, (1-68) 

and 



7esc(EQ ^ T c 2 

5 ^ E? ' 



« d-69) 

Th 



where the escape rate is given by 



7e 



(1.70) 



7 2 21n(2coshx)' 
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Then we find from Eqs. (1.59,1. 62): 

Wi„ = T c 2 ff (J±-^ . (1.71) 
One can see that due to the fact that both the cooling rate W ou t given by 




U, a.u. 



Figure 1. 7. A sketch of the Coulomb blockade peak shape in the dynamic localization regime 
without taking into account the phonon cooling (solid line): at small U < t/ m in the dephasing is 
dominated by the electron escape (peak), at larger U - by electron-electron collisions (plateau), 
and finally, at U > Umax the cooling is insufficient, the dynamic localization is destroyed, and 
the dot is in the Ohmic regime. The Ohmic curve is also shown for reference by the dashed line. 

Eq.(1.67) and the energy pumping rate W in given by Eq.(1.71) are proportional 
to T c 2 ff , the effective electron temperature drops out of the global energy bal- 
ance Wi n = W ut- As the result, the dc conductance of the dot is remarkably 
independent of the gate voltage: 
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Note that since G/Gq < 1 such a plateau solution exists only when 



At the same time for large enough detuning U > U m i n ~ T* the condition 
Eq.(1.69) is fulfilled because the escape rate 7 CS c(^) decreases and the ef- 
fective temperature increases with increasing U. However, with U and T e ff 
further increasing the DL regime Eq.(1.68) is destroyed and G{U)/Gq starts 



G(U) dependence is shown in Fig.7. 

5. Summary 

In conclusion, we have considered the manifestation of dynamic localization 
in a solid-state system - a quantum dot under time-dependent perturbation. We 
have shown that in a closed dot without electron interaction the dynamic local- 
ization is possible and we developed the theory of weak dynamic localization 
for an arbitrary time-dependent perturbation. Then we considered a realistic 
case of a dot with electron-electron interaction weakly connected to leads. We 
have demonstrated that if the dot is subject to a periodic perturbation the dy- 
namic localization can result in a plateau at the tail of the Coulomb blockade 
peak in the linear dc conductance vs the gate voltage and we established the 
conditions under which such a plateau arises. 
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